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Various models for pion multiplicity distributions produced in relativistic heavy ion colli¬ 
sions are discussed. The models include a relativistic hydrodynamic model, a thermodynamic 
description, an emitting source pion laser model, and a description which generates a nega¬ 
tive binomial description. The approach developed can be used to discuss other cases which 
will be mentioned. The pion probability distributions for these various cases are compared. 
Comparison of the pion laser model and Bose-Einstein condensation in a laser trap and with 
the thermal model are made. The thermal model and hydrodynamic model are also used 
to illustrate why the number of pions never diverges and why the Bose-Einstein correction 
effects are relatively small. The pion emission strength g of a Poisson emitter and a critical 
density r] c are connected in a thermal model by r]/n c = e~ m ' T < 1, and this fact reduces 
any Bose-Einstein correction effects in the number and number fluctuation of pions. Fluctu¬ 
ations can be much larger than Poisson in the pion laser model and for a negative binomial 
description. The clan representation of the negative binomial distribution due to Van Hove 
and Giovannini is discussed using the present description. Applications to CERN/NA44 and 
CERN/NA49 data are discussed in terms of the relativistic hydrodynamic model. 

PACS numbers: 24.10.Jv, 21.65.+f, 24.85.+p, 25.75.-q 


I. INTRODUCTION 

Pion multiplicity distributions produced in heavy ion col¬ 
lisions are of current interest for several reasons. First, 
pions are by far the main component of the produced par¬ 
ticles in very high energy heavy ion collisions. For exam¬ 
ple, in the SPS experiments at CERN hundreds of pions 
are produced, and this number may go up considerably at 
RHIC energies. In the Landau hydrodynamic model, the 
number of pions scales at S' 1 / 4 where S = E 2 m . Secondly, 
pions are a useful tool for studying HBT effects coming 
from Bose-Einstein correlations. HBT two-particle cor¬ 
relation experiments give information about source pa¬ 
rameters of the emitting system. Thirdly, if many pi¬ 
ons are produced by a strongly emitting source at high 
enough density, a pion laser may be formed according 
to one model owing to Bose-Einstein symmetrization ef¬ 
fects. Pions with a nonzero chemical potential can also 
show the phenomena of Bose-Einstein condensation at 
some critical density. Since Bose-Einstein condensation 
of a finite number of atoms in a harmonic oscillator or 
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laser trap has recently been seen, the possibility of seeing 
effects associated with the formation of a pion laser or an 
enhanced condensation into a ground state would be in¬ 
teresting. A comparison of Bose-Einstein condensation 
of atoms in a harmonic oscillator trap and the pion laser 
model will be given in this paper. Fourthly, the possi¬ 
bility of intermittency behavior in pion distributions has 
been a concern for over a decade. A distribution that is 
widely used to discuss intermittency is the negative bino¬ 
mial distribution, and a model which generates a negative 
binomial distribution will be developed in the framework 
of the present description. The clan representation for 
the negative binomial distribution due to Van Hove and 
Giovannini will also be discussed. This representation is 
a useful way of characterizing experimental data. 

In this paper we present a unified description of several 
currently discussed models for evaluating pion distribu¬ 
tions. Our main focus will be on four particular cases, but 
other possibilities exist and will be discussed briefly. The 
examples we will consider in detail are 1) a thermody¬ 
namic description, 2) a relativistic hydrodynamic model, 
3) an emitting source model, and 4) a model which gives a 
negative binomial distribution. These cases cover a wide 
spectrum of possibilities. Thermodynamic descriptions 
have been extensively used [|j (s| for the production of 

pions in heavy ion collisions and also are discussed in 
the quark matter meetings [|]j|] . Emitting source models 
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are reviewed in Refs. |],|7j. The emitting source model 
that we discuss in detail is originally due to Pratt || and 
discussed further in Refs. Application of the neg¬ 

ative binomial distribution to high energy collision data 
can be found in Refs. (TT) (T^j, and its connection with 
intermittency are discussed in Refs. Hydrody¬ 

namic descriptions are also frequently applied to heavy 
ion collisions, and their descriptions appear in the quark 
matter proceedings, Refs. |[|||. The relativistic hydro- 
dynamic model gives a good description of single and 
double inclusive data from Cern experiments carried out 
by the CERN/NA44 and CERN/NA49 collaborations. 
We will use an approach that initially was developed for 
cluster distributions JItJ Q, but applications to cy¬ 
cle class problems were note d [|17| and further studied 
in more detail in Refs, [^lj Cycle class problems 

appear in Bose-Einstein and Fermi-Dirac statistics be¬ 
cause the density matrix must be symmetrized or an¬ 
tisymmetrized in Feynman’s path integral approach to 
statistical mechanics [|d). A detailed application of the 
approach discussed below to Bose-Einstein condensation 
phenomena can be found in Refs. || 22 |j' 2 .'i| |. The connection 
between cluster models and cycle class problems is shown 
in Fig. 1. Specifically, a cycle of length k corresponds to 
a cluster of size k, and the number of cycles of length k 
corresponds to the number of clusters of size k. 
Application of results from cluster models to the theory of 
disordered systems can be found in Refs. Some of 

the results of this reference were connected to probability 
distributions which appear in Quantum Optics. Methods 
from Quantum Optics have been used by Weiner j27j] 
p9| et al. to discuss particle production and correlations 
in heavy ion collisions. 


II. PION PROBABILITY DISTRIBUTIONS 

In this section we will discuss several models of pion pro¬ 
duction in heavy ion collisions. These models, which are a 
thermodynamic description, a relativistic hydrodynamic 
model, an emitting source model, and a negative bino¬ 
mial model, can be described using a unified formalism. 
Before giving specific results, we will give some general 
expressions that we will use. These expressions were ini¬ 
tially developed for the description of the fragmentation 
of nuclei into clusters. However, because of the corre¬ 
spondence between cycles and clusters shown in Fig. 1, 
the cluster results can be used to describe cycle class 
problems, and the pion production models that we con¬ 
sider can be re-expressed in such terms. This correspon¬ 
dence has already been illustrated in the description of 
Bose-Einstein condensation. 

This section is divided into various subsections. First, 
we give some general results. Then we will apply these 
results to the cases mentioned above. Comparisons will 


be made for the various cases and with other situations 
already discussed such as Bose-Einstein condensation. 


A. General Results 


In describing a fragmentation process, partitioning prob¬ 
lems, or cycle class problems, a weight is given to each 
division of A = Yk^ nk - Fig ure 2 illustrates the possi¬ 
ble divisions of A = 4. A vector n = (ni, ri 2 ,..., n k , ua) 
specifies the division where nk is the number of clusters 
with k nucleons, cycles of length fc, etc. A general block 
picture for n is shown in Fig. 1, and the specific block 
pictures of A = 4 are shown in Fig. 2 along with cluster 
and cycle class divisions and harmonic oscillator excita¬ 
tions. The type of weight IF 4 (n.x) given to any ft that 
we will consider can be written as 


W A (n, x) 


1 

Z A (x) 


n 


X 


nfc 

k 


n k \ 


(1) 


Here nQ- are Gibbs factorials, Z A (x) is the canonical par¬ 
tition function, and x k is a function of thermodynamic 
quantities such as volume F, temperature T. This x k was 
called a tuning parameter in Ref. and a specific form 
Xk = x/k was studied in detail in Refs. |l 7[|l s| . Other 
forms for x k were used to study percolative features of 
simple statistical model [|2l|j and Bose-Einstein conden¬ 
sation of atoms in a box | 22 f| and atoms in a laser trap 
[22| . The form of x k for the pion models to be discussed 
will be given below. 

Given WA(n, x) of the type of Eq. (1), various quantities 
are easily calculated. For example the ensemble averages 
of n k , given by (n k ) = Y W A (n, x)n k , where the sum is 
over all partition of A, is simply |l7lj2G|| 


(n k ) = x k 


Z A ~k{x) 

Za{x) 


( 2 ) 


This result can be easily understood heuristically by not¬ 
ing that the n k times the Gibbs factorial 1/rifc! gives 
1 /{nk — 1)!. This operation removes one column of k 
blocks from Fig. 1. The sum over all partitions of A, 
which we call tta(ji) in Y Wa{ n, x)n k , then gives rise to 
the factor Z A - k {x). The x k factor in Eq. ( 2 ) is from 
the shift x^ k = x k x r ^ k ^ 1 where the x^. k ~ l is the appro¬ 
priate power of Xk for the partitions with one column of 
k blocks removed. The Z A (x) in Eq. (2) is present in 
Wa(H, x) and is in the final answer. Since the constraint 
A = Y k k nk always applies to any partition of A, then it 
is also true when written as A = Yk k i nk )- Substituting 
the result of Eq. (2) into this constraint equations gives 
A = Yk kxkZA— k / ZA^ and rearranging terms results in 
the recursive relation 


Za{x) 


1 

A 


kx k Z A -k(x) 

k 


( 3 ) 
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with Zq = 1. This recursive equation is a very simple 
way of obtaining Za(x), and it applies to any Wa{^,x ) 
of the form of Eq. (1). To use it, the Xk& have to be 
specified. As an example Xk = x/k gives 0 .® 

Z A (x) =x(x + 1)...(x + A- 1 )/A! = r ^ A f (4) 

The Za(x) is the canonical partition function of a sys¬ 
tem of A objects ( A nucleons, A atoms, ...). The 
grand canonical partition function z(x,u) is obtained 
from the canonical partition function Za(x) or the weight 
Wa(ti,x). For example, the generating function of 
W A {x) = E IA( ii) W A (n,x) is @ 


z(x, u) = exp [ux i + u 2 x 2 + u 3 x 3 + ...] = Wa(x)u a 

A=0 

( 5 ) 

The u is taken to be u = in thermodynamics, where 
p is the chemical potential. For /i = 0, u = 1 and 


z(x) = exp 


H xk 

.k =1 


( 6 ) 


The ratio 


Za(x)/z(x) = Pa{x) 


( 7 ) 


In a thermal model or in a description based on statisti¬ 
cal mechanics, the distributions associated with pions are 
obtained from the assumption that an equilibrium is es¬ 
tablished in some interaction volume V at some temper¬ 
ature T. Specifically, statistical thermodynamics gives 
the Bose-Einstein occupancy factor 

/(«> = ( 10 ) 

for an energy level at e and (3 = 1/ksT. This result 
is based on the grand canonical ensemble with the pion 
chemical potential p = 0. For pions in a box of volume 
V, the thermal properties of pions can be obtained from 
the density of states factor Vd 3 p/h 3 and from the energy 
e, momentum p , relation e = (p 2 + m 2 ) 1 / 2 . In particular, 
with h = c = ks = 1 , the following relations can be 
obtained for thermal pions: 



k= 1 


l = ((E/V)+p)/T 


is the ratio of the canonical partition function to the 
grand canonical partition function. The mean and vari¬ 
ance of A follow once the Xk’s are specified: 

(A) = Y kx k ( 8 ) 

k 

and 

(A 2 ) - (A ) 2 = Y, k2 *k ■ (9) 

k 

These last two equations follow from Eq. (5) by differen¬ 
tiating this equation with respect to u and setting u = 1 . 
These expressions form the main general results to be 
used in this paper. Other results can be found in Refs. 

0-HI- 


B. Pion Distributions in a Thermal Model 

We begin this subsection by first giving some well-known 
results on the thermal properties of an ideal gas of rela¬ 
tivistic pions. Some of these results will then be used to 
motivate some ideas, developments, and comparisons to 
be made later. Moreover, since the thermal description 
is extensively used, many readers will recognize these re¬ 
sults and this is helpful for the subsequent development 
of the choice of Xk for this particular model and other 
models to follow. 


Here N/V is the total number density, E/V —the energy 
density, pV —the equation of state, and S/V —the en¬ 
tropy density of pions. The various K 's are MacDonald 
functions. The ultrarelativistic photon-like limit m —> 0 
orT>m and the nonrelativistic limit e = p 2 /2m+m can 
be obtained from these equations using the following lim¬ 
its for the MacDonald functions A'„ that appear. With 
z = km/T : K v {z) = (T(^)/2) j (1/2 z) v for 2 : —> 0 and 
K v (z) = yj{ k/Zz) e~ z for z —> 00 . Here T(i/) is a gamma 
function. In the above Eq. (11), the sum over k gives the 
corrections arising from Bose-Einstein or BE statistics, 
with the k = 1 term the Maxwell Boltzmann or MB limit. 
Specifically, the k = 2,3, ...00 give the enhancements 
coming from BE statistics. For example, for N/V the 
ratio N/V/{N/V)mb = C(3) in the massless pion limit. 
Here (1V/V)mb is the k = 1 term only. The £(3) = 1.2 is 
the Riemann zeta function. The results of Eq. (11) arise 
when the occupancy factor of Eq. (10) is expanded as a 
power series /(e) = e~^ e (l + e~^ e + e -2/3e + ...). The 
integrals over (Vd 3 p/h 3 )e~ k/3e with e = \Jp 2 + m 2 give 
rise to the various A'’s. 

The probability of having nk pions in level fc, given that 
the mean number in that level is ( nk ) = 1/ (e^ £fe — l) is 
the geometric distribution [|o| 


p(n k ) = 


(n fc p 

«n fc ) + lp +1 


( 12 ) 
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The (n k ) = ^ p{n k )n k and (n k ) - ( n k ) 2 is 
n/c= 0 

(nl) - (n k ) 2 = (rife) (1 + (■ n k )) (13) 

The Poisson limits is (n k ) — (n k ) 2 = ( n k ) and is real¬ 
ized when MB statistics applies. The above expression 
on fluctuations pertain to a particular energy level. We 
now turn to an evaluation of the fluctuations in the to¬ 
tal number of pions and its associated probability dis¬ 
tribution. Such quantities are important in an event- 
by-event analysis of data. The probability of A-pions 
is given by the ratio of the canonical partition function, 
Zn, to the grand canonical partition function z, with 

OO 

z = Zn and Zo = 1. To obtain Zn,z we use results 
TV—0 

from Sec. II.A with u = = 1. For a thermal model 

the choice of x k is 

VT 3 /ton 2 1 /ton 

Xk= 2^P) *** a (*r) (14> 

Later, we will give other choice of x k - That this is the 
correct choice follows from (A) = kx k where (A) /V 
in the grand canonical ensemble is given by Eq. (11). 

Also 2 = e PA T = eP kXk where PV/T = J2 x k follows 

k 

also from Eq. (11), and z = e EkXk was a result developed 
in Sec. II.A. The result for 2 = e^ Xk can be expanded to 
give the Wn(x) of Eq. (5) and its associated Wjv(n, x) = 
]”[ x^ k /n k \/Z n(x) of Eq. (1). The Wn (n, x) is the weight 
given to n = (?ii, ri 2 ■ • • ua) in a cycle class decomposition 
of a permutation. For example the permutation 


fl 2 3 4 5 6 7\ 
^1324756/ 


(15) 


has two cycles of length 1 which are 1 —> 1 and 4 —> 4, 
one cycle of length 2 which is 2 —> 3 —» 2, and one cycle of 
length 3 which is 5—>7—> 6 —>5. The cycles appear in 
Feynman’s path integral density matrix approach in Sta¬ 
tistical Mechanics p4| , and permutations appear because 
the density matrix must be symmetrized. An example of 
this approach to the ideal Bose gas in a box [J22| and 
to the ideal Bose gas in a harmonic oscillator trap |2^ j 
shows the utility of the recurrence relation of Eq. (3) in 
obtaining the canonical partition function Zn ( x ) and all 
its associated thermodynamic quantities. We will briefly 
summarize these two cases in the next section. 

Using this choice of x k . the probability of A pions is 
P N = Zas/z where Zn is obtained from Eq. (3) with 
x k given by Eq. (14). The fluctuations (A -2 ) — (A ) 2 = 
Y^k 2 x k . An alternative method for obtaining Pn was 
recently given by Becattini, Giovannini, and Lupia [pl| . 
They used projection methods to obtain Pn which in¬ 
volves taking A-derivatives of z with respect to u and 


setting u = 0 . Specifically 


Pn = 


1 


jN 


lim 


N\ w-*o dw N 


exp 


V 

2i r 2 


J d 3 plog^l-e £/,T wj 


(16) 

Table 1 summarizes Pn for T = 163 MeV, V = 22.5, and 
m/T = 0.85. This choice generates the same Pn distri¬ 
bution given in Ref. J3^|. Here the Pat’s are calculated by 
the simple recurrence relation of Eq. (3). Computation¬ 
ally, Eq. (3) is very easy to work with since computers 
evaluate recursive relations quickly as well as the sums 
for (A) and (A 2 ) — (A ) 2 that appear in Eq. ( 8 ). 

Two limiting behavior of x k are of interest. In the ultra- 
relativistic limit or massless pion limit, the x k of Eq. (14) 
is simply 


V , 1 

** = ¥ 

and in the nonrelativistic limit 

T , — ^ r ~km/T 1 

Xk ~\ 3 T G kP 2 


(17) 


(18) 


where A t = h/(2n mT) 1 / 2 . The factor e~ km P is the 
Boltzmann factor cost function for creating fc-pions, and 
this factor arises from the mass of the pion. We will 
introduce a quantity t, where x k ~ l/k T . In cluster 
models r is Fisher’s critical exponent |32j. In Eq. (17) 
r = 4 and in Eq. (18) r = 5/2. 


C. Relativistic Hydrodynamic Model and 
Applications to CERN/NA44 and CERN/NA49 
Experiments 

We also calculated the pion probability distribution us¬ 
ing a relativistic hydrodynamic model: HYLANDER-C 
[|28| . This particular model belongs to the class of mod¬ 
els which apply 3+1-dimensional relativistic one-fluid- 
dynamics. It provides fully three-dimensional solutions 
of the hydrodynamical relativistic Euler-equations [ 351. 
HYLANDER-C and its earlier version, HYLANDER [ 331, 
have been successfully applied to various heavy-ion reac¬ 
tions at SPS energies. HYLANDER-C and HYLANDER 
were used to reproduce J3(| simultaneously mesonic and 
baryonic rapidity and transverse momentum spectra of 
the S + S reaction at 200 AGeV. Corresponding mea¬ 
surements have been performed by the NA35 Collabo¬ 
ration |uj]. Based on the successful description of the 
measured single-inclusive spectra, predictions for Bose- 
Einstein correlation (BEC) functions j3f]|F| were made. 
Those predictions turned out to agree quantitatively with 
the measurements [ [Io|j4T| . The model also reproduces the 
photon data for S+Au collisions at SPS energies jl2| and 
gives a simple explanation for the “soft-px puzzle” 
and the complex behaviour of the radii extracted from 
pion and kaon correlations and explains the difference in 
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the extracted radii for pions and kaons in terms of a cloud 
of pions due to the decay of resonances which surrounds 
the fireball (pion halo) [f38f . 

In the following, we discuss some further results |34j] 
for 158A GeV Pb + Pb collisions measured by the 


CERN/NA44 144 § and CERN/NA49 Collaborations 


J46|. The by HYLANDER-C reproduced Pb+Pb data 
have been obtained while using an equation of state 
(EOS) with a phase transition to a quark-gluon plasma 
at a critical temperature Tc = 200 MeV ( cf. Refs. [ fl7| , 
p9| and Refs, therein). This EOS does not depend on 
the baryon density, and thus the freeze-out energy den¬ 
sity translates directly into a fixed freeze-out temperature 
Tf. The choice for the freeze-out temperature is in these 
calculations, Tf = 139MeU. 

The number of thermal pions ( 7 r + , and 7r°) that are 
produced in the model calculations is 853. In particu¬ 
lar, this number does not include any resonance decay 
contributions. But this total number of thermal pions 
includes the Bose-Einstein enhancement. To be specific, 
the thermal pion spectra are calculated by evaluating the 
integral 


dN 1 f 

p^da^x/) 


d 3 p (2tt)3 

|VX(4)1 

- 1 

exp 

[ TfK) \ 


(19) 


where r = d+ 1. Then 0 = expfaj^r)], (IV) = XdC(d) 
and (TV 2 ) - (TV) 2 = x d ((d - 1) with £(n) = J2kLi I/*"- 
Some values of £(ra) that often appeal in statistical me¬ 
chanics are £( 1 ) = 00 , £( 2 ) = 7 t 2 /6 = 1.645, C(3) = 1 . 202 , 
((4) = 7t 2 /90 = 1.082, £(3/2) = 2.612. In two dimen¬ 
sions (TV 2 )-(TV) 2 = 00 so that the fluctuations of mass¬ 
less particles becomes large. 

For nonrelativistic atoms in a box p^| in d-dimensions 


x 1 


( 21 ) 


with r = d/2 + 1, x = L d / Ay, A t = h/{2'KinT) 1 / 2 . Con¬ 
densation occurs in d = 3 at x = x c = A/£(3/2) = V/\j. c 
in the limit A, V —> 00 . Manifestation of this Bose- 
Einstein condensation are 1) a power law distribution of 
cycle lengths plj.pTl : (n k ) ~ 1/fc 5 / 2 , 2) the specific heat 
Cy has a cusp |30|, and 3) the ground state occupany 
nc.s. becomes macroscopically large Q: nc.s./A = 
1 - x/x c = 1 - (T/T c ) 3 / 2 , for T < T c . The T c is de¬ 
termined by the density A/V from A/£(3/2) = V/X^. 
For d = 1, 2 T c —> 0 so that there is no sudden conden¬ 
sation at a nonzero T. The occupancy of any level with 
energy e k is 


n k 


E 


—0Ckk ZA—k 

Z A 


( 22 ) 


where da M is the differential volume element and the in¬ 
tegration is performed over the freeze-out hypersurface, 
E. u IJ (x) and Tf are the four velocity of the fluid element 
at point x and the freeze-out temperature, respectively. 
p 11 = ( E,p ) is the four-momentum of a pion. 

By expanding the Bose-Einstein occupation factor in 
terms of a Maxwell-Boltzmann factor and a series of 
higher order corrections, we can obtain the distribution 
of (Efc’s which are summarized in Table 2. 

Using these results for x k , the fluctuation in the ther¬ 
mal pions is (TV 2 ) — (TV) 2 = ]T) k 2 x k = 312.97 or a 20% 
enhancement over the Poisson result of Xi = 260. The 
probability distribution associated with the thermal pi¬ 
ons is obtained from P/v = Z n /z where Zn is given 
by Eq. (3) and 2 is given by Eq. ( 6 ). Figure 3 shows 
the result for P N , which is a perfect Gaussian distribu¬ 
tion (solid line) with mean (TV) = 283.9 and variance 
(TV 2 ) - (TV) 2 = 312.97. 


For A —> 00 , ZA-k/ZA —^ e /3/ifc and thus we obtain the 
well-known occupancy factor 


1 

Uk ~ — 1 


(23) 


As T —> T c , p, —> eo = 0. 

For atoms in a laser trap given by a harmonic oscillator 
well with frequency u> and level spacing hu; [Q 


_ 1 x kd ' 2 
k k (1 — x k ) d 


(24) 


where x = e For small x 


_ x kd ! 2 (k B T\ 
Xk k d+1 \ huj ) 


There is no condensation at nonzero T for <2=1, but 
there is condensation for d = 2, 3 at T c 7 ^ 0. 


D. Some Special Examples 

In this subsection, we will illustrate some of the general 
results of Sec. II.A with three specific examples. For 
massless particles in a box of sides L in d-dimensions 
(neglecting spin), the x k is 

L d (k B T\ d 1 _ x d 
k 2 d ~ 3 TT d ~ 1 y he J k T k T 


E. Emitting Source, Pion Laser Model 

A simple and exactly solvable emitting source models for 
pions was first introduced by S. Pratt [||. The model 
has an interesting feature. Namely, a Poisson emitter of 
pions can behave like a pion laser when Bose-Einstein 
symmetrization effects are included. A manifestation of 
this behavior is a large enhancement in the number of 
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pions. This enhancement is into a zero momentum state 
and is similar to Bose-Einstein condensation of atoms in 
a box which also condense into the zero momentum state. 
An important result of this condensation is a reduction 
in the intercept in the two particle correlation functions. 
Such effects were discussed in Refs. ]27j Q. Recently 
Csorgo and Zimany [ |l0| have given simple exact expres¬ 
sions for Pratt’s model. Some simple results were also 
presented by Chao, Gao, and Zhang Q. This subsection 
presents some further discussions of this model using re¬ 
sults from Sec. II. A when applied to this case and shows 
some similarities with the Bose-Einstein condensation of 
atoms in a laser trap as given by Eq. (24). 

First, we will give expressions for Xfc. Once Xk is given, 
then all the general results of Sec. II.B apply. In Pratt’s 
model, 77 is the mean pion multiplicity generated by a 
Poisson emitter. Then, in the Poisson limit, the prob¬ 
ability of having TV pions is Pjv = {jl N / AH) e~ v . This 
is also the result of the thermal model of Sec. II.B 
in the Maxwell-Boltzmann limit, where 77 = x\ = 
(VT 3 /2tt 2 ) (m/T)K 2 {m/T). Corrections to this result 
obtained from Bose-Einstein symmetrization effects can 
be developed by evaluating a quantity Xk (called Ck in 
Ref. 0). The Ck are referred to as combinants ju|. The 
importance of combinants can be found in the work of 
Gyulassy and Kaufman |l|,[40 . While the Xk of the ther¬ 
mal model is given by Eq. (14), in the emitting source 
model Xk is given by evaluating 


xk= T f n n 

i =1 i—1 

X II s (2 + ft ) ’ Xi ) ( 25 ) 

i=1 ' ' 

where po = P = Pk and x\ = 77 . We neglect any temporal 
dependence to keep the discussion simple. 

The s(p, x) is the source strength. As in Ref. [|| this 
source strength is taken to be a thermal Gaussian of the 
form 


S(p,x) 


_£_ e -P 2 /2mT+x 2 /2R 2 

{{2tt) 2 R 2 mT) 3/2 


(26) 


For this choice of s(p, x), the evaluation of the x^’s can 
be done analytically 0. We now give a discussion of the 
liming behavior of Xk before giving the general behavior. 
For small fc and R 2 /2 l/8mT (typically R 2 /2 ~ 6 and 
l/ 8 mT~ 1/6) 


Xk 


( R 2 mT) 


3/2 


(. R 2 mT ) 3/2 



(27) 


Thus Xk ~ 1/fc 4 and r = 4 for small fc. For large fc 



(28) 


Now Xk ~ 1/fc and r = 1. The p c is given by 

/ __ 3 / 2 

? 7 c = ( R 2 mT + VR 2 mT + — j (29) 


Table 3 gives some values of Xk for R 2 /2 = 6 and 
l/ 8 ?uT=l/ 6 . 

The general behavior for Xk was explicitly written by 
Csorgo and Zimany ( 10 , and is 


Xk 


f 1 

k ( k/2 k/2\ 3 

£ 7 + - 7 - ) 


(30) 


with 7 ± = 1/2 (l + y ± </l -(- 77 ) and y = 2R 2 mT — 1/2. 
The 77 c = in this notation. We rewrite this as 



to illustrate its close relation to atoms in a laser trap as 
given by Eq. (24). For d = 3, Eq. (24) gives 


x 3k /2 

%]$ — q" 

fc (1 — x k ) 


(32) 


with x = e - nu l kBT . 

Once the Xfc’s are given, the z = exp[J]xfe], ( N) = 
Y^kxk, (N 2 ) — {N) 2 = Ylh 2 Xk and Pn = Zn/z follow 
with Zjy given by the recurrence expression of Eq. (3). 
Since Xk = l/k(rj/ri c ) k for large fc, (N) —> 00 when 
77 > 77 c . At 77 = 77 c the (N) ~ ]T) 1/fc which diverges loga¬ 
rithmically. We now compare these results with thermal 
models. In thermal models, Xk is given by Eq. (14) and its 
nonrelativistic and ultrarelativistic limits in Eq. (17) and 
Eq. (18). No divergence in (TV) occurs in either the ultra¬ 
relativistic limit where (N) = (V/i r 2 ) T 3 (( 3), or the non¬ 
relativistic limit where (TV) = (K/Ai/) < 73/2 ( m/T ) with 

S3/2(w/T) = J/(e _m / T ) fe /fc 3 / 2 . Note, however, that 
fluctuations, (TV 2 ) — (TV)" ', in the two-dimensional ultra¬ 
relativistic case discussed in Sec. II.C can be infinite. To 
further illustrate the difference, we note that in the ther¬ 
mal model all x^’s are given by Eq. (14) so that xi = 77 is 
not independent of V, T, m where V ~ R 3 . Table 4 shows 
values of Xk for V = R 3 with R 2 /2 = 6 and mT = 3/4. 

Using 77 = xi, the Xk can be rewritten as Xk ~ p k and the 
results are also shown when written in this form. The re¬ 
sults of Table 4 are very similar to Table 3. However, the 
large Xk limit is (V/ (A^ • fc 5 / 2 )) e ~ krn / T , If we take x\ to 
be this result with fc = 1, the Xk = (l/fc 5 / 2 ) {j]/p c ) k with 
7 ?c = U/Ay ~ (R 2 mT ) 3 ^ 2 , which is similar to Eq. (29). 
Note that 77 = r] c e~ m ^ T so that r]/r] c = e -m / T < 1. Con¬ 
sequently, the sum x^fc = (TV) always converges. Thus 
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there can be no result (N) —> oo. The Boltzmann factor 
e -m/T SU pp r esses the property (N) —> oo. This hidden 
connection between 77, 7 ?c, and r\lr] c < 1 is not present 
in the pion laser model which treats 77 as independent of 
?7 C or X\ as independent of the other Xk s. The thermal 
model has this connection and shows why no divergence 
in ( N) ever occurs in the thermal model. Similar conclu¬ 
sions also apply to the relativistic hydrodynamic model. 


F. Negative Binomial Distribution 

As we have seen, pion probability distributions are de¬ 
termined once the choice of the Xk s is specified. The 
thermal equilibrium model specifies a particular form for 
Xk as does the emitting source model. Bose atoms in a 
box or laser trap also can be discussed in terms of Xk- 
In this subsection we will present a choice for Xk which 
results in a negative binomial distribution. The choice 
given below is partly motivated by results from the other 
models. The negative binomial is a probability distribu¬ 
tion that is frequently used to characterize pion yields 
p| pi . It has also been used in the discussion of in- 
termittency p|,^6|. Van Hove and Giovannini JTTJ m 
has given a clan model representation of it which will also 
be mentioned. 

The various forms for Xk that we have encountered so far 
in this paper and also in the application of this approach 
to cluster yields suggest looking at a behavior for Xk given 

by 

x k =ay k /k T (33) 

For example, the nonrelativistic thermal model has a = 
V/A|., y = e -m / T , and r = 5/2. The ultrarelativistic 
thermal model has a = yyT 3 , y = 1, and r = 4. The 

emitting source model, for small k, is a = ( R 2 tttT) 3 / 2 , 
y = r]/(R 2 mT ) 3 / 2 , r = 4 and is a = 1, y = ?7/?7c, t = 1 
for large k. We will now show that the choice 

Xk = ay k /k , (34) 


Y[(y k ) nk = y N since fen*, = N , where the Y\{y k ) nk 
k k 

appears in Wa (n. x) of Eq. (1). One way to obtain this 
result for Zn is directly from the recurrence relation of 
Eq. (3). The mean number of pions can be obtained from 
P kxk and is 

(N) = -W- (35) 

1 -y 

The variance follows from (-V 2 ) — (TV) 2 = ay /(l — y) 2 or 
(N 2 ) - (N) 2 = (N) (l + ^\ (36) 


The result of Eq. (36) is the result for the variance of a 
negative binomial distribution. This result is also similar 
to that of Eq. (13) when a = 1, but Eq. (13) refers to 
fluctuations in a particular level k, while Eq. (36) has 
no reference to a level. The fluctuations can be much 
larger than Poissonian because of the factor (1 + ( N ) /a). 
To see that the probability distribution associated with 
this model (xk = ay k /fc) is indeed a negative binomial 
distribution, the Pn = Zn/z can be easily evaluated 
and rewritten as 


P N = 


N + a-1 
N 



( ( N)/a \ 


(37) 


The result of Eq. (37) is the standard form of the negative 
binomial distribution. 

An interesting transformation of the negative binomial is 
to a set of clan variables N c , n c , given by pj[ 


N c = a log ( 1 + 
N r . 


(N) 


(38) 


The N c equals the average number of clans, and n c is the 
average number of particles in a clan. The variables N c , 
n c are a useful way of analyzing experimental data. In 
terms of y and z = 1/(1 — y) a 


which has r = 1 , generates a negative binomial prob¬ 
ability distribution. This way of generating a negative 
binomial distribution can be found in Ref. [ 5C ], more re¬ 
cently in Ref. |H], and was used in Refs. [17] and [ po| . 
The a and y will be related to the mean number of pions 
and to their fluctuations. We will also connect them to 
the clan variables of Van Hove. Using results from Ref. 
0 and Sec. II.A, the following behavior for various rel¬ 
evant quantities are obtained. First, the grand canonical 
partition function z = exp (J/ Xk) = 1/(1 — y) a which can 
easily be verified by noting ay k /k = —a log(l— y). The 

k 

canonical partition function Z n can be obtained and is 
p| Zn = y N T(a+N)/T(a)N \—see also Eq. (4). The ex¬ 
tra y N that appears here is from the y k in Xk ■ Note that 


N c = log 2 

n c = —7 -- . ^ . -r (39) 

(1 — 2/) log(l — 2 /) 

Thus n c only depends on y and N c is simply connected 
to z. Note that in thermodynamics (see also Sec. II.B), 
the grand canonical partition function z is connected to 
the equation of state. Specifically, pV/ksT = log ,3 so 
that pV = N c ksT if this connection is made. 


G. Other Models 

Besides the choice r = 1, other choices for r can be 
solved. Many of these other choices were studied in the 
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framework of cluster models [0 , but the results obtained 
also apply to particle production of bosons. Some discus¬ 
sion of these issues can be found in Ref. |2{| where results 
from a somewhat general theory of disordered systems 
were connected to distributions used in particle physics 
and quantum optics. 

III. CONCLUSIONS 

This paper discusses several models of pion production in 
a unified way by showing that they are different choices 
of a quantity which we called Xk ■ The models considered 
are a thermal description of pion production, a relativis¬ 
tic hydrodynamic model, an emitting source pion laser 
model, and a model that gives rise to a negative binomial 
distribution. The results of these models are compared 
with each other. Comparisons are also made with Bose- 
Einstein condensation of Bose particles or atoms trapped 
in a harmonic oscillator well and in a box of d-dimensions. 
The quantity Xk appears in the weight given to each cycle 
class decomposition of the symmetric group. The index 
k is the length of a cycle. The cycle classes appear when 
the density matrix is symmetrized as in Feynman’s path 
integral approach to statistical mechanics. The quan¬ 
tities Xk, k = 1,2,3,... completely determine the be¬ 
havior of the system as discussed in Sec. II.A. Specifi¬ 
cally the mean number of pions and their fluctuations are 
just moments of the Xk distribution. The grand canoni¬ 
cal partition function z is an also simply obtained from 
these Xk s, and the canonical partition function Zjq can 
be obtained by a simple recursive procedure which con¬ 
tains the Xk- The probability of ./V-pions is Pjv = Zn / z. 
These methods and relations were initially used by one 
of us in cluster problems 0 n and subsequently ap¬ 
plied to Bose-Einstein cycle class problems 0 , 0 . They 
were also used in a discussion of the pion laser model 
in Refs. § 0- The two problems, cycle class and 

clusters, can be mapped into one another as shown in 
Fig. 1. Specifically the index k is either the cluster size 
or cycle length, and nk is either the number of cycles of 
length k or the number of clusters with k nucleons. After 
presenting various general expressions, specific cases are 
discussed and compared. We started with the thermal 
model because of its long history and its widespread use 
in heavy ion collisions.The thermal model was used to 
motivate some of the quantities, such as the x^s that 
also appear in the other cases that are considered. The 
role of the Boltzmann factor e~ km / T in Xk, which rep¬ 
resents the cost function for making fc-pions, is impor¬ 
tant in the thermal model. The mean number of pions, 
their fluctuations, and the probability distribution for N- 
pions are discussed and compared with Poisson results. 
The results for a relativistic hydrodynamic model are also 
discussed. The relativistic hydrodynamic model is used 
in the analysis of CERN/NA44 and CERN/NA49 data. 


The probability distribution for thermal pions is shown 
to be a Gaussian distribution. The fluctuations in pi¬ 
ons is about 20% above the Poisson limit. This result 
of slightly enhanced fluctuations is also consistent with 
thermal model results. 

The emitting source model of S. Pratt is next discussed. 
It has an interesting feature of the possibility of large 
departures from Poissonian results. This model also has 
some other interesting properties that are not present in 
the thermal and hydrodynamic models. Namely, a crit¬ 
ical pion density r] c exists and, if the Poisson emission 
strength r) equals this density, an infinite number of pi¬ 
ons results from Bose-Einstein enhancement factors. The 
thermal and hydrodynamic models lack this property be¬ 
cause rj ~ r\ c e~ rn l T . Some result of the emitting source 
model are shown to be similar to Bose-Einstein conden¬ 
sation in a laser trap when this trap is taken to be a 
three-dimensional harmonic oscillator well. Another dis¬ 
tribution which can have large fluctuations above Pois¬ 
son results is the negative binomial distribution. The 
clan representation of the negative binomial distribution, 
introduced by Van Hove and Giovannini, is discussed in 
terms of quantities that appear in the specification of the 
Xk s. The clan representation is a useful way of charac¬ 
terizing data, and it is compared with the thermal model 
results. 

This work has been supported by the U.S. Department 
of Energy. 
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TABLE I. The probability Pn of A pions compared to a 
Poisson distribution with the same (A). Also given are the 
cycle class parameters Xk obtained from Eq. (14). 

A P N 

Poisson 

k 

Xk 

0 .311 

.293 



1 .346 

.359 

1 

1.10812 

2 .208 

.221 

2 

.04812 

3 .0903 

.0904 

3 

.00588 

4 .0320 

.0278 

4 

.00107 

5 .00998 

.00682 

5 

.00024 

6 .00287 

.00139 

6 

.00006 


TABLE II. The distribution of Xk in a 
dynamic Model 

Relativistic Hydro- 

k 

Xk 

k 

Xk 

1 

260 

6 

.0067 

2 

9.957 

7 

.0016 

3 

1.047 

8 

.0004 

4 

.163 

9 

.0001 

5 

.031 




TABLE III. Values of Xk- Unless rj is large, the Xk fall very 
rapidly and thus the {N) and (A 2 ^ — (A) 2 are governed by 
the first few k' s. For the coice R 2 /2 = 6 and 1/8 mT = 1/6, 
the r/c = 42.875. However, the second Xk = X 2 is comparable 
to at r] = 432, the X3 to X2 at 77 = 140, the xa to X3 at 
rj = 90, etc. 


k 

Xk 

k 

Xk 

1 

V 

6 

? 772.43 x 10 1U 

2 

if / 432 

7 

rf /1.32 x 10 12 

3 

?7 3 /6.07 x 10 4 

8 

?7 8 /7.40 x 10 13 

4 

?7 4 /5.47 x 10® 

9 

7773.8 x 10 15 

5 

?7 5 /3.91 x 10 s 

10 

?7 lo /2.07 x 10 17 
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TABLE IV. Values of Xk in a thermal model using the same 
choice of R 2 , T as in Table 2. The Xk is also written as rf/bk 
where xi = r) and bk is a number that is given in the table. 
The results rf /bk are very similar to that of Table 3. 


k 

Xk 

Xk = V k /bk 

1 

20 

V 

2 

0.991 

?7 2 /404 

3 

0.142 

rf/ 6.21 x 10 4 

4 

0.309 x 10" 1 

r? 4 /5.17 x 10 6 

5 

00.835 x 10~ 2 

rf/ 3.83 x 10 8 

6 

00.257 x 10 -2 

rf/ 2.57 x 10 10 

7 

000.869 x 10~ 3 

rf/1A7 x 10 12 

8 

000.312 x 10~ 3 

? 7 8 /8.20 x 10 13 

9 

000.118 x 10~ 3 

77 9 /4.34 x 10 15 

10 

0000.462 x 10' 4 

rf°/2.22 x 10 17 


FIG. 1. Cycles and Clusters. In the mapping of a cycle 
class problem into a cluster problem and vice versa, the cycle 
length corresponds to the cluster size; the number of cycles of 
length k corresponds to the number of clusters with k nucle¬ 
ons. Also shown is the vector n = (m, ri 2 ,..., ua) as a block 
diagram. The number of vertical blocks is k , and the number 
of columns with k blocks is nk . The total number of blocks 
is A. 


FIG. 2. Various Parallels. Various parallels are shown 
which include the partitioning of an integer, a correspond¬ 
ing block picture, a cluster, a cycle class representation of 
each partition, and an equivalent harmonic oscillator excita¬ 
tion. Some places where these various columns appear are as 
follows: The partitions of an integer appears in counting the 
number of irreducible representation of the symmetric group. 
Block pictures and rotated versions of the ones shown ap¬ 
pear in combinatorial analysis (Ferrer’s block diagram) and 
in group theory (Young tableaus). Cluster problems appear 
in many areas of physics such as nuclear multifragmentation, 
percolation, randomly broken objects, and the fragmentation 
of atomic clusters. They also appear in group structure such 
as in the distribution of city sizes, etc. Cycle class problems 
appear in Bose-Einstein and Fermi-Dirac statistics, speckle 
patterns in phase space, Feynman’s theory of the A-transition 
in liquid He, in random permutations, etc. Excitations in 
the harmonic oscillator appear when evaluating the micro- 
canonical partition function of Bose atoms in a laser trap. 
Another application is to the Veneziano-Hagedorn mass spec¬ 
trum, where each excitation is an elementary particle. 


FIG. 3. Results for Pjv for thermal pions, n~, in 158 
A ■ GeV Pb+Pb collisions. The solid line uses contributions 
k —> 600 (already k = 10 gives the fully converged result). 
The dotted lines only uses up to k = 2 terms. In the figure, 
we have (solid line) E Xk = 271.2, (N) = E kxk = 283.9, and 
(N 2 ) - (N) 2 = E k 2 x k = 312.97. 
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